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Abstract 


The  performance  of  on  airborne  laser  terrain  mapper,  which  collects 
slant  range  and  reflectance  data  simultaneously,  is  determined  in  part 
by  the  spatial  filtering  effects  of  the  finite  laser  beam  spot  size. 

An  analysis  of  the  spatial  filtering  effects  :uni  the  resulting  distortion 
of  the  terrain  profile  slant  range  to  tue  aircraft  is  made  to  interpret 
the  interaction  between  the  terrain  characteristics  and  the  resulting 
laser  line  scan  information.  Lmphasis  is  placed  on  the  range  interactions 
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subcarrier  phase  modulated  by  the  terrain  heigh:  profile.  "his  phase 
modulated  signal  is  low  nass  filtered  to  represent  the  effects  of  the 
spatial  filter.  Numerical  analysis  is  utilises  to  evaluate  this  nonlinear 
filtering  problem  ana  to  examine  the  range  error  in  terms  of  the  maximum 
range  height  variations.  The  resulting  mean-squared  range  distortion 
error  is  added  to  the  range  squares  distortion  error  originating  from  the 
receiver  noise  in  the  Phase  Locked  Loop  phase  tracking  circuit.  The 
distortion  sue  to  spatial  filtering  is  shown  to  increase  with  the  laser 
modulation  frequency  while  the  error.;  u.e  to  detection  noise  decrease 
with  frequency.  Thus,  the  performance  of  the  receiver  can  be  optimized 
according  to  the  terrain  maximum  range  height  variations  ana  the  desired 
resolution  by  designing  tiie  properly  sized  optics  and  choosing  an  optimum 
frequency  of  modulation.  In  addition,  the  general  spatial  filtering 
model  is  extended  to  evaluate  the  performance  of  a forward  looking  laser 
terrain  mapper  using  Low  grazing  angles. 
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SPATIAL  FILTERING  DESIGN 


CONSIDERATIONS  FOR  A LASER  LINE 

SCANNING  SENSOR 

I . Int roduction 

The  design  of  an  airborne  terrain  mapping  laser  scanner  is  being 
studied  by  the  Air  Force.  This  proposed  sensor  would  collect  range  ana 
reflectivity  data  which  could  be  combined  and  used  to  construct  three 
dimensional  imagery  of  the  terrain.  Some  of  the  possible  applications 
are  referenced  in  an  AGARD  report  (Ref  l).  Figure  1 is  a block  diagram 
of  the  laser  line  scanning  sensor  which  consists  of  the  following  basic 
components:  laser  illuminator,  line  scanner,  direct  detection  receiver, 
optics  train,  signal  processing  electronics,  display  device,  and  a 
recording  device.  The  line  scanner  is  used  to  scan  the  terrain  in  small 
strips  perpendicular  to  the  aircraft's  forward  motion,  while  the  line 
scans  are  incremented  along  the  y-axis  by  the  aircraft's  forward  motion 
as  depicted  in  Figure  2.  The  beam  intensity  pattern  as  represented 
along  the  x-axis  is  b(x),  and  the  width  of  this  intensity  pattern  along 
the  x-axis  is  the  beam  spot  size  which  is  designated  by  X.  The  laser  is 
intensity  modulated  by  a periodic  signal.  The  laser  is  used  to  illu- 
minate the  terrain  with  this  periodic  intensity  modulation  and  this 
signal  is  detected  by  a direct  detection  receiver.  Within  the  receiver, 
a particular  harmonic  of  the  periodic  signal  is  phase  tracked,  and  the 
resulting  phase  is  compared  to  the  phase  of  the  transmitted  signal.  The 
phase  difference  between  the  two  signals  is  converted  to  the  slant  range 
from  the  aircraft  to  the  terrain  spot  being  scanned.  Since  the  phase 
used  to  determine  the  slant  range  is  relative  phase. 


it  may  contain 


Mod  2tt  ambiguities.  Thus,  the  range  measured  is  relative  range  and  will 
contain  a c/2fm  range  ambiguity,  where  c is  the  velocity  of  light  and 
is  the  frequency  of  the  detected  harmonic. 

Problem 

The  purpose  of  this  thesis  is  to  analyze  the  spatial  filtering 
effects  of  the  finite  beam  size,  and  to  determine  the  interaction 
between  the  various  laser  scanning  sensor  parameters  and  the  terrain 
characteristics.  This  analysis  of  the  spatial  filtering  can  be  used  to 
determine  the  sensor's  fundamental  limits  of  resolution  The  primary 
goal  of  this  thesis  is  to  examine  the  image  degradation  caused  by  the 
spatial  filtering.  In  addition,  there  is  an  interest  in  extending  the 
analysis  of  the  spatial  filtering  to  the  analysis  of  the  spatial  filtering 
for  the  design  of  a forward  looking  sensor  utilizing  small  grazing  angles. 
Previous  work  done  by  Robinson  (Ref  2)  included  an  analysis  of  the  sensor 
with  sinusoidal  scanner  modulation  which  determined  the  conditions  for  no 
distortion  of  the  terrain  input  signals  by  the  spatial  filtering.  These 
results  form  the  main  approach  to  the  problem  and  will  be  reviewed  in  the 
next  chapter.  Also,  research  has  been  conducted  on  several  promising 
types  of  temporal  modulation  for  the  laser  (Ref  3),  but  the  analysis 
needs  to  incorporate  the  effects  of  spatial  filtering  to  be  complete. 

Such  an  overall  analysis  will  be  used  to  optimize  the  sensor  performance 
and  will  impact  future  sensor  designs. 

Assumptions 

In  order  to  limit  the  scope  of  the  problem,  several  assumptions  are 
made  in  modeling  the  spatial  filtering  for  the  laser  line  scanning  sensor. 
First,  the  receiver  is  a direct  detection  receiver.  Second,  only  periodic 
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temporal  intensity  modulation  or  modulation  which  ran  be  represented  as 
a Fourier  Series  will  be  considered.  Third,  the  beam  will  be  modeled 
as  having  only  one  dimension  which  will  be  along  the  scanner  direction 
of  motion  (x-axis  in  Figure  l).  In  turn,  the  terrain  wilL  only  be  modeled 
along  the  same  axis.  Fourth,  the  receiver  is  aligned  and  functioning 
properly.  Fifth,  the  receiver  noise  is  wideband  noise.  Finally,  it  is 
assumed  that  the  post  detection  signal  processing  will  not  further 
distort  the  desired  information. 

Thesis  Overview 

The  analysis  of  the  spatial  filtering  problem  starts  with  a review 
of  the  spatial  filtering  model.  A general  approach  for  modeling  the 
terrain  as  a statistically  stationary  process  is  examined.  i'he  terrain 
range  and  reflectivity  are  modeled  as  Gaussian  Random  Processes.  Pue 
to  the  nonlinear  nature  of  the  input  r.-uige  signal,  only  general  relation- 
ships between  the  terrain  character i st. ics  and  the  sensor  paramet.ors  for 
no  distortion  are  derived.  Because  the  general  problem  i s beyond  tin- 
scope  of  this  study , some  special  cases  are  examined  ami  analysed. 
Specifically,  the  cases  which  are  .uiaLyv.ed  are  for  constant  terrain 
range  ,-uid  for  linearization  of  the  phase  modulated  range  input  signal. 

One  other  special  case  which  is  examined  is  the  case  for  sampling, 
receivers.  This  analysis  applies  to  pulse  modulat ion  systems. 

The  following  chapter  is  composed  of  an  analysis  for  spatial 
filtering  distortion  for  a single  sinusoidal  range  input  ah  ng  with 
constant  reflectivity.  fho  distortion  due  to  the  spatial  filtering  is 
modeled  as  an  ideal  low  pass  filter.  I’he  range  error  caused  by  the 
distortion  from  the  spatial  filtering  is  added  to  the  range  error 
caused  by  tin  distortion  from,  the  receiver  noise  in  the  phase  tracker  ef 
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the  receiver  giving  the  total  range  error.  The  resulting  distortion  data 
is  plotted  in  terms  of  the  terrain  range  characteristics,  sensor  beam- 
width,  frequency  of  modulation,  and  the  range  error.  The  distortion 
results  are  applied  to  compute  the  range  error  for  several  specific 
cases  including  an  extension  to  the  forward  looking  case. 

In  the  remaining  chapter  the  general  results  are  summarized  and 
some  recommendations  for  solving  the  general  spatial  filtering  problem 
for  any  type  of  range  input  are  presented. 
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I I . Spatial  Filter  Model 

A one-dimensional  spatial  filtering  model  representing  the  beam 
filtering  process  of  the  terrain  range  and  reflectivity  of  the  laser 
line  scanner  was  developed  by  Robinson  and  Reinmtui  (Ref  J,  li ) . This 
model  applies  only  for  periodic  sinusoidal  intensity  modulation  of  the 
laser,  but  some  crude  generalizations  can  be  made  for  other  types  of 
periodic  intensity  modulation.  The  model  was  originally  used  to  deter- 
mine the  conditions  under  which  the  terrain  range  and  reflectivity  signals 
would  not  be  distorted.  A short  review  of  the  model  along  with  the 
conditions  obtained  for  distortionless  received  signals  are  presented  in 
the  following  sections.  This  leads  to  the  ['rime  empnasis  of  tnis  tiiesis 
which  is  to  analyze  the  distortion  of  the  range  ;uia  reflectivity  signals 
caused  by  the  beam  filtering.  Because  the  general  case  of  spatial 
filtering  is  not  easily  analyzed,  it  is  necessary  to  model  and  idealize 
the  description  of  the  terrain  range  find  reflectivity,  and  the  be  sun 
filtering  function. 


The  received  signal  at  the  detector  output  is  represented  by  the 
following  integral  (Ref 


s ( t ) = P T p(x)  h(x-vt)  f(t  - ^r^)dx  (i) 

where  r(x)  is  the  slant  range  from  the  terrain  to  the  aircraft,  b(x)  Is 

the  beam  filtering  function,  p(x)  is  the  terrain  reflectivity,  v is  the 

scan  velocity,  f(t)  is  the  temporal  intensity  modulation,  c is  the 

velocity  of  light , and  i'  is  the  peak  optical  power  which  includes  the 
* o 

propagation  and  optics  losses . Ttiis  one  dimensional  representation  of 
the  spatial  filtering  is  valid  if  the  two  dimensions l beam  filtering 
function,  b(x,  y),  is  extremely  narrow  along  the  v-nxis,  or  if  the  range 
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r(x)  is  treated  as  an  average  of  the  range  values  in  the  y dimension 
over  the  beam  spot  size.  Eq  (l)  is  valid  for  all  types  of  intensity 
modulation.  It  is  convenient  to  assume  that  f(t)  is  periodic  and  has  the 


following  Fourier  Series  representation: 


£ c „e 


+J  lnnt't 


m“,  where  f i. 

m 


the  fundamental  frequency  and  the  c^'s  are  the  Fourier  Series  coeffi- 
cients. By  inserting  the  Fourier  Series  representation  of  f(t)  in  Eq  (l), 
it  can  be  rewritten  as  follows: 


s(t)  = P0  _£  cnr  p ( x ) b(x-vt)  e+J‘-7rnlm(  t_  0 (Ref 

Q CO  il  — uO 


) (2) 


The  Fourier  Transform  of  Eq  (2)  can  be  determined  using  two  Fourier 
Transform  theorems,  and  is  the  expression  in  Eq  (3), 


3(f)  = PQ  n=L  ITT  D(I^ilS) 


(Kef  2:3)  ( 3 ) 


where  B(f)  is  equal  to  the  Fourier  Transform  of  b(x),  and  D(f)  is  equal 

to  the  Fourier  Transform  of  p(x)  exp  ( J^TTn.!m.r.^ x. ) . For  harmonics  with 

c 

nonzero  c^'s,  the  periodic  modulation  shifts  the  spectra  so  that  is  is 
centered  at  multiples  of  the  fundamental  frequency.  The  amount  of  infor- 
mation received  is  independent  of  the  scan  velocity  because  both  D(f) 
and  B(f)  are  scaled  by  v.  It  is  important  to  realize  that  the  bandwidth 
of  one  spectrum  alone  does  scale  with  v as  this  will  impact  the  post, 
detector  electronics  design.  Although  all  the  spectral  contain  range  and 
reflectivity  information,  the  information  will  vary  between  harmonics 
because  the  range  ambiguity  interval  is  different.  The  constants  in 
front  of  Eq  (3)  will  be  ignored  because  they  do  not  contribute  to  the 
range  or  reflectivity  data.  This  Leads  directly  to  the  spatial  filter 
model  which  Robinson  derivea  (Ref  2)  and  which  is  shown  in  Figure  3. 

The  model  in  Figure  3 is  the  equivalent  complex  model  for  the  detection 
of  any  harmonic,  n f ^ . One  important  consideration  is  that  the  range 
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data  being  measured  is  relative  range  and  may  contain  Mod  2n  ambiguities 
depending  upon  A , the  fundamental  wavelength  of  tne  temporal  modula- 
tion; and  the  range  variations  witnin  the  beam.  Unless  otherwise 
stated,  this  model  will  be  used  throughout  the  remainder  of  the  thesis. 
The  first  approach  to  the  spatial  filtering  analysis  is  the  stocastic 
modeling  of  the  terrain  parameters. 

Gauss ian  Random  Processes  Model 

An  original  approach  to  the  spatial  filtering  problem  is  to  model 
the  terrain  range  and  reflectivity  as  Gaussian  Hanuom  Processes  tGRP), 
compute  the  power  spectral  density  of  the  input,  ana  determine  how  tne 
Fourier  Spectrum  is  distorted  by  the  spatial  filter.  The  following 
assumptions  are  used  to  model  the  terrain  as  a GRP: 

1.  p(x)  and  r(x)  are  statistically  independent  GRP  with  r.ero  rae.uis. 

2.  The  terrain  scene  is  statistically  stationary  so  tiiat  the  auto- 
correlation functions  R (x,x')%  the  range  autocorrelation  function,  arid 
Rp(x,x'),  the  reflectivity  autocorrelation  function,  can  be  written  as 
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R (d)  and  R (d)  respectively,  where  d = x-x' . 
r p 

By  using  the  above  assumptions  and  the  characteristic  function,  it  can 
easily  be  shown  that  the  autocorrelation  function  of  p(x)exp(-i-J^,I-nIr::1  _) 


R (d)  = R0(d)exp[-  (l!^£si)2(R  (0)  - R ( d) ) ] (U) 

OP  err 

The  Fourier  Transform  of  Eq  (k)  is  the  power  spectral  density  and  can  be 

written  as  follows : 

/Unnfm,2Rr(0)  ■ Uirnfm,  2Rr(d) 

c:  ( ' z ' f Id  r ' c ' 1 Cel 


Sd(f)  = e"v  — 


Fd[Rr(d)]»  Fd[e 


where  Fd  is  the  Fourier  Transform  operator  on  the  variable  d.  If  the 
bandwidth  of  the  Fourier  Transform  of  Rp(d)  is  3p , and  the  bandwidth  of 

the  Fourier  Transform  of  exp(  ( t'r^^j  qs  3 5 then  crudely  the 

c * 

bandwidth  of  S,(f)  is  the  sum  of  these  two  banawidtns.  Thus  for  the 
d 

input  signal  to  be  filtered  with  no  distortion, 

bp  + V < Bf  (6) 

<P 

where  3^  is  the  spatial  bandwidth  of  the  beam  filter.  Depending  upon 
the  characteristics  of  the  range  and  reflectivity  autocorrelation  func- 
tions, it  is  usually  very  difficult  to  compute  3 and  B . Another 
bandwidth  measure  is  mean  squared  bandwidth  of  the  power  spectral  density 
(Ref  5:100).  It  can  be  used  to  rewrite  Eq  (o)  as  follows: 


[ W~  + < [bT 

nI  P >1  r^  >1  t (7) 

where  is  the  reflectivity  mean  squared  bandwidth,  B^  is  the  mean 
squared  bandwidth  of  the  Fourier  Transform  of  exp(  ( Br^^),  and 

is  the  mean  squared  bandwidth  of  the  filter.  Eq  (7)  can  be  rewritten 


in  terms  of  the  beam  spot  site,  the  mean-squared  bandwidth  of  the  range 


i 


I 


< 

I 


process,  B2,  using  some  previously  derived  results  (Ref  2:5-o): 


. B2 

>J  •> 


1 ^ ^ 7t  n 1 

u r"1 


) R ( o ) B2 
r r 


1 

X2 


(8) 


In  h'q  (8)  R^(0)  is  the  mean  squared  range  variation,  and  l/X2  (X  is  the 

beam  spot  sine)  is  approximately  equal  to  the  mean  squared  bandwidth  of 

~2  ~ 

the  spatial  filter.  B is  generally  much  larger  than  B . For  no  dis- 

‘ $ r 

tortion  and  the  input  bandwidth  Bt  being  small,  can  be  proportionally 

large  compared  to  the  case  for  B^  being  large.  Similarly,  if  the  input 

. 2 . i ■» 

bandwidth  Br  is  small.  Bp  can  be  large  compared  to  the  case  for  Bj  being 
large.  Thus,  there  is  an  interaction  between  r(x)  and  p(x)  which  depends 
on  the  respective  input  bandwidths  and  directly  affects  the  distortion  of 
the  resulting  output  signals. 


Distortion  Effects 

The  GRP  model  is  good  for  determining  the  conditions  for  no  distor- 
tion, but  it  is  substantially  troublesome  to  use  for  determining  how  the 
signal  will  be  distorted  if  the  terrain  input  signal,  p ( x ) exp ( '4-iTn-— 1 ■— , 

has  a bandwidth  which  exceeds  the  bandwidth  of  the  spatial  filter. 
Typically,  this  condition  is  encountered  for  rough  terrain  that  has  sharp 
edges,  and  rapidly  varying  terrain  heights  and  changes  in  reflectivity, 
and  will  result  in  some  distortion  of  the  terrain  range  and  reflectivity 
signals.  The  distortion  process  caused  by  the  beam  filtering  is  a linear 
process;  however,  the  phase  modulated  range  input  signal  is  nonlinear, 
thus  tlie  resultant  output  signals  which  are  filtered  will  be  nonlinear. 
This  problem  is  the  standard  amplitude-phase  modulated  signal  which  Is 
filtered;  and  is  studied  by  Middleton  (Ref  o:oj8-oi2),  and  many  others 
in  the  communications  field.  Unfortunately,  this  problem  does  not  yield 
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any  exact,  tractable  analytic  solutions,  thus  most  solutions  are  generally 
arrived  at  by  numerical  means.  However,  there  are  some  specialized 
results  which  are  made  possible  by  modeling  the  terrain  range  variation, 
the  terrain  reflectivity,  and  the  beam  function  for  idealized  cases. 

Linear  Approach 

One  approximation  which  considerably  simplifies  the  spatial 
filtering  problem  is  the  series  expansion  of  the  exponential  and 
truncation  of  this  series  to  include  just  the  first  two  terms.  How,  the 
terrain  range  input  to  the  filter  is  linear  in  r(x)  and  the  spatial 
filtering  is  equivalent  to  the  low  pass  filtering  of  the  product  of  two 
amplitude  modulation  processes.  The  linear  model  representation  of  the 
beam  spatial  filtering  is  shown  in  Figure  4.  For  the  linear  approxima- 


s(x) 

■“  LA/V 

Fig.  4.  Linear  Approximation  for  Spatial  Filtering 

4imfrar(x) 

tior.  to  be  valid,  the  following  inequality  must  hold,  - <<  1. 

Typically,  this  means  that  the  left  side  of  the  inequality  must  be  less 
than  it/12,  or  the  phase  variations  are  less  than  15°.  The  linear  case 
ensures  that  the  bandwidth  of  the  phase  modulated  process  is  narrowband. 
If  the  variations  in  range  are  small,  then  for  no  distortion  the  spatial 
frequency  variations  of  the  range  can  be  considerably  larger  than  when 
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compared  to  the  general  phase  modulation  case.  When  the  n=0  term  is 
detected,  there  is  a reflectivity  term  without  any  dependence  on  r(x). 

If  this  reflectivity  term  could  be  detected,  it  could  be  used  to  remove 
the  reflectivity  dependence  of  the  other  terms,  and  if  done  properly 
could  decrease  the  spatial  bandwidth  on  the  input  signal  to  the  spatial 
filter.  (This  discussion  on  the  removal  of  the  reflectivity  dependence 
is  applicable  to  the  general  problem,  not  just  the  linear  approximation . ) 
Constant  Range 

Another  method  for  analyzing  the  beam  interaction  with  the  terrain 
reflectivity  is  to  assume  that  the  terrain  range  is  constant,  so  that 
Eq  (l)  can  be  written  as: 

s(t)  = T° b(x-vt)  f ( t— c i ) p(x)ax  (9) 

— oo  X 

where  c^  is  a constant  range  term.  The  Fourier  Transform  of  Eq  (9)  is 
the  convolution  expression  in  Eq  (10),  where  F(f)  is  the  Fourier  Transform 

S(f)  = e-jClF(f)#  ~ B*( f/v)R( f/v ) (10) 

of  f(t),  3(f)  is  the  Fourier  Transform  of  b(x),  and  R(f)  is  the  Fourier 
Transform  of  p(x).  The  whole  process  can  be  modeled  in  the  time  domain 
as  shown  in  Figure  5*  Therefore,  for  the  constant  range  assumption,  onJ.y 


X =vt 


Fig.  5-  Constant  Range  Spatial  Filter  Model 
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the  reflectivity  variation  will  be  observed.  Also,  the  distortion  which 
will  be  observed  is  simply  the  low  pass  filtering  of  an  amplitude 
modulated  signal.  This  implies  that  the  sensor  will  be  able  to  resolve 
reflectivity  variations  along  the  x-uxis  down  to  the  sice  of  the  beam 
spot  size.  The  entire  constant  range  discussion  is  analogous  to 
observing  the  Fourier  Series  component  for  n=0  as  previously  discussed 
under  the  linear  approximation. 

Another  Represent  at  Ion-Sampling 

The  effects  of  spatial  filtering  on  the  terrain  signals  can  be 
characterized  for  sampling  receivers.  This  is  particularly  valuable 
because  the  sampling  concepts  can  be  used  to  extend  the  spatial  filtering 
results  from  the  case  of  sinusoid  modulation  to  pulse  types  of  periodic 
modulation  to  pulse  types  of  periodic  modulation  for  the  scanner. 
Initally,  the  terrain  range  will  be  assumed  to  be  constant,  so  that  only 
the  variations  in  reflectivity  will  be  observed.  Using  the  constant 
range  representation  which  was  previously  presented,  the  beam  spatial 
filtering  can  be  represented  as  shown  in  Figure  6,  where  p'(x)  is  the 
filtered  reflectivity.  For  a pulse  modulated  system,  p'(x)  would  be 


Fig.  6.  Reflectivity  Spatial  Filtering  for  Constant  Range 
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sampled  at  the  pulse  repetition  rate,  and  the  sampled  p'(x)  would  be 
used  to  reconstruct  p'(x).  To  reconstruct  p'(x),  it  is  necessary  to 
determine  how  close  p'(x)  should  be  sampled  along  x,  or  what  the  pulse 
repetition  rate  should  be.  Since  p’(x)  is  bandlimited  by  b(x)  which 
has  an  approximate  bandwidth  of  l/X,  then  p'(x)  shoula  be  sampled  at  2/X 
or  twice  every  beamwidth  to  prevent  aliasing.  If  this  type  of  argument 
is  expanded  to  include  the  terrain  range  input,  the  detected  signal  will 
also  be  bandlimited  by  the  spatial  filter  to  l/X. 

This  representation  can  be  extended  to  look  at  the  more  general 
problem  which  encompasses  other  types  of  periodic  modulation  for  the 
scanner.  This  problem  can  be  examined  by  taking  the  inverse  Fourier 
Transform  of  Eq  (3),  which  can  be  represented  by  a sum  of  signals  which 
are  demodulated  to  baseband  and  passed  through  sampling  filters  and  is 

given  by  Eq  (ll).  In  principle,  Eq  (ll)  could  be  applied 
co  x A*r(x)nfm 

s(x)  = Z c p ( x ) e J c #b(x)  (ll) 

n 

to  determine  the  output  signal  by  using  an  infinite  number  of  receivers, 
where  each  receiver  detects  a different  spatially  filtered  harmonic. 

From  this  modeling  of  the  sampling  receiver  two  conclusions  can  be 
drawn.  First,  for  a constant  range  sampling  receiver,  the  reflectivity 
variations  which  are  on  the  order  of  the  spot  size  will  be  resolved 
providing  that  the  spatial  filter  output  is  sampled  at  2/X.  Similarly, 
for  the  sampling  receiver  which  simultaneously  detects  reflectivity  and 
range  variations,  the  minimum  detectable  range  or  reflectivity  varia- 
' tions  along  the  x-axis  will  be  several  beamwidths  depending  upon  the 
bandwidth  of  the  input  signal  to  the  spatial  filter.  Clearly,  the 
distortion  of  the  input  terrain  signal  is  a function  of  its  bandwidth. 

I 
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deootid , no  natter  what  pulse  repetition  rate 


chosen,  assuming  that 


III.  Distortion  Analysis  and  Results 


Within  this  chapter,  the  distortion  which  results  from  the  spatial 
filtering  is  discussed  ana  analysed.  As  previously  mentioncu,  the  dis- 
tortion is  caused  by  the  bandlimiting  of  the  spatial  filter.  The  distor- 
tion process  is  of  considerable  importance  in  analyzing  the  laser  line 
scanning  sensor  because  the  resolution  of  the  sensor  is  limited  by  the 
distortion  process.  Unless  a measure  of  distortion  is  refined,  the  term 
distortion  is  rather  nebulous.  For  the  purpose  of  this  thesis,  distortion 
is  defined  in  terms  of  the  mean  squared  error  of  the  detected  signals. 

For  other  measures  of  distortion  and  their  calculation,  tne  appropriate 
communication  literature  can  be  consulted  (Ref  7-15).  After  a measure  of 
distortion  is  chosen,  somehow  the  distortion  must  be  computed.  Unfor- 
tunately, the  modulation  process  is  a nonlinear  process,  and  thus  the 
distortion  can  only  be  computed  by  numerical  analysis  with  the  aid  of  a 
digital  computer.  Once  the  distortion  has  been  modeled,  it  is  desirable 
to  compute  the  distortion  in  terms  of  the  sensor  and  the  terrain 
parameters.  After  this  step  is  accomplished,  the  distortion  uata  can  be 
used  to  optimize  the  sensor  parameters  for  minimum  distortion,  given  a 
particular  set  of  terrain  parameters. 

Ueneral  Analysis  of  the  .Spatial  Filtering  Distortion 

To  make  a general  analysis  of  the  spatial  filtering  distortion,  it 
is  desirable  to  express  the  distortion  in  terms  of  an  input-output 
relation.  Although  this  representation  simplified  the  problem,  it  does 
not  provide  an  analytical  solution  for  modeling  the  distortion  effects. 
Figure  7 is  the  illustration  of  the  spatial  filter  representation  for 
the  distortion  input-output  relation.  Each  output  signal,  p'(x)  and 
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Specialized  Distortion  Analysis 

The  general  problem  can  be  solved,  but  the  solution  requires  an 
extensive  numerical  analysis,  and  the  solutions  are  very  specialized. 

Therefore,  some  simplifications  and  idealizations  are  necessary  to  make 
the  problem  more  tractable  and  to  provide  some  useful  results.  Probably, 
the  simplest  idealization  is  to  assume  that  the  linear  approximation  is 
valid.  The  outgrowth  of  this  condition  is  that  the  phase  modulated 
range  signal  is  narrowband  and  this  portion  of  the  signal  behaves  as  an 
amplitude  modulated  signal.  Thereby,  the  linear  case  is  analogous  to 
having  an  input  signal  which  is  the  product  of  two  amplitude  modulated 
signals.  The  advantage  of  the  linear  case  is  that  the  distortion  process 
for  amplitude  modulation  is  analytically  tractable  and  well  understood. 

The  linear  case  can  be  implemented  for  the  sensor  by  sufficiently 
decreasing  the  fundamental  frequency  of  modulation  until  the  phase  varia- 
tion is  less  than  15°.  Unfortunately,  when  the  fundamental  frequency  is 
decreased,  the  accuracy  of  the  relative  range  decreases  also.  Another 
consequence  of  the  linear  approximation  which  also  applies  to  the 
general  problem  is  that  at  the  center  frequency  of  the  output,  where  n=0, 

| 

the  reflectivity  term  is  decoupled  from  the  range  term.  This  reflectivity 
term  could  be  separately  detected  and  used  to  separate  the  reflectivity 
dependence  from  the  input  or  output  of  the  spatial  filter. 

Two  methods  are  presented  for  removing  the  reflectance  dependence. 

One  method  decouples  the  reflectivity  at  the  input,  while  the  second  scheme 
decouples  the  reflectivity  term  at  tne  output. 

1.  Decoupling  the  reflectivity  term  at  the  input  to  the  spatial 
filter  can  be  implemented  by  using  a two  beam  scanner  where  the  first  beam 
illuminates  the  terrain  with  a cw  signal  and  is  used  to  recover  the  f 
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reflectivity  signal.  The  reciprocal  of  the  reflectivity  signal  is  used 
to  amplitude  modulate  the  second  beam  which  follows  directly  behind  the 
first  beam.  Thus,  the  reflectivity  dependence  is  nearly  removed  :uid 
only  the  phase  modulated  range  signal  is  filtered.  This  method  is  expen- 
sive to  implement  because  of  the  additional  equipment  required,  but  the 
output  signal  is  completely  decoupled  from  the  reflectivity. 

2.  Decoupling  the  reflectivity  at  the  output  can  be  accomplished 
by  using  a d.c.  receiver  pickoff  and  dividing  the  dependency  of  this 
signal  out  of  the  existing  signal.  However,  there  is  one  serious  consid- 
eration, the  distortion  terms  will  still  be  present  even  though  the 
reflectivity  dependence  has  been  removed.  These  distortion  terms  could 
cause  errors  in  the  detected  terrain  range  signal. 

There  are  two  other  specialised  cases  which  are  considerably  simpler 
to  analyse,  constant  range  and  constant  reflectivity.  The  constant  range 
case  is  of  little  practical  consequence;  but  it  does  illustrate  a partic- 
ular case,  namely,  the  lowpass  filtering  of  :ui  :unplitude  modulated  signal 
where  the  distortion  process  is  well  understood,  in  the  case  of  constant 
reflectivity,  it  is  now  reasonable  to  accomplish  a numerical  analysis  of 
the  distortion  because  the  input  terrain  signal  is  strictly  a chase 
modulated  signal,  and  the  resulting  :uialysis  is  much  easier.  Using  tills 
assumption  of  constant  reflectivity,  the  range  distortion  wiil  be  modeled 
and  computed. 

Range  Distortion  Mode  Ling, 

The  distortion  of  the  range  signal  can  be  easily  modeled,  and 
numerical  results  can  be  obtained  by  assuming  that  the  reflectivity  is 
constant  and  that  the  terrain  range  can  be  modeled  as  a periodic  function. 
The  modeling  of  the  terrain  reflectivity  as  a constant  is  valid  for 


slowly  varying  reflectivity  or  when  the  reflectivity  dependence  is 
removed  as  previously  illustrated.  The  periodic  function  representa- 
tion of  r(x)  is  an  appropriate  comparison  to  the  real  terrain  over  short 
distances  along  the  x-axis  which  could  include  up  to  several  beamwidths. 
Actual  calculations  of  the  distortion  to  the  range  signal  will  include 
additional  assumptions  in  order  to  simplify  the  numerical  analysis. 

Fourier  Series  Representation  of  the  Phc.se  Modulated  Periodic  Range 

Function 

A periodic  function  which  is  phase  modulated  can  be  represented  as 
a product  of  Bessel  Fourier  Series.  The  terrain  range  which  is  assumed 
to  be  a periodic  function  has  a Fourier  Series  representation. 


oo  + j t !*• 

ooale  > where  the  a^'s  are  the  complex  Fourier  coefficients. 


and  f is  the  fundamental  spatial  frequency  of  the  range  function.  The 

K 

phase  modulated  range  function  can  be  expressed  as  follows  using  the 

trigonometric  representation  of  a complex  exponential: 

Unnfnirtx)  Wf^J^a  exp[+J2  fk*x] 

e J c = e c c 

A .Unnfmr(x)  ^ ^ , 

e ‘ c = „Tr  exp[+j  t1Tn  -n-  a«(cos  2nf,.Jlx  + jsin  2nf  .4x)]  (13b) 

)C“ — oo  c iv  iv 


(13a) 


Using  the  Bessel  function  identities  of  a complex  exponential  with  a sine 
or  cosine  argument,  the  right  side  of  Eq  (13b)  can  be  rewritten  as  an 
infinite  product  of  Bessel  Series  as  shown  in  Eq  (lE),  where  n^  equals 


a*!e 


j.  j ^ffnfmr(  x) 


Ermf 


n Z (1  + Jll)J  (^ 

= _OOii=_od  ° n r» 


4=_a>,[=_ 


)e 


+jq(2irf..lx  + 04) 


(1U) 


For  exajnple,  if  the  riinge  is  a single  sinusoid,  where  r(x)  = r^sintJnf  x)  , 
then  its  Bessel  Fourier  Series  representation  is: 
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+jUffnfmrmsin(2l,fr*) 
e o 


) ,.2irfrxq 

= t J,(^Hlrm)e+J 

q=-“  4 c 


If  (■—  1 m)  is  replaced  by  tf  (which  is  usually  interpreted 

iation  index),  L'q  (15)  can  be  rewritten  as  a real,  plus  an 
infinite  series  of  the  following  form: 


(15) 


as  the  rnodu- 
i magi nary 


+J  tJsinPnf  ,.x 
e 1 


= Jo(tf)  + 2[Jo(b)cos2(r(2fr)x  + J1+(ti)cos  2n(Ufr)x 
+ 2j  [ Jj_(  3)sin27rfrx  + J.(b)sin  2n(if  )x  + 


This  form  of  the  phase  modulated  range  signal  can  be  easily  used  to 
evaluate,  by  means  of  numerical  analysis,  the  spatial  filtering  distor- 
tion . 

Sensor  Distort  ion  Mode  1 i ng  for  a o i nusoi  u Range  Ir.j  ut 

For  ease  of  computation,  the  terrain  range  signal  is  assumed  to  be 
a single  frequency  periodic  sinusoid  of  the  form  r(x)  = r sin( jnf  x) , 
where  r is  the  peak  amplitude  range  variation  and  f^  is  the  spatial 
frequency  of  the  terrain  range  input.  Two  sources  of  distortion  are 
considered  in  the  modeling  of  the  range  error  for  the  laser  line  scanning 
sensor.  Besides  distortion  being  caused  by  the  spatial  filtering,  there 
will  be  some  range  error  induced  by  the  receiver  noise  within  the  Phase 
Locked  Loop  (PLL)  circuit  which  is  used  to  track  the  phase  of  tne  detected 
signal.  The  two  distortion  processes  are  independent  and  each  interacts 
with  the  range  signal  as  illustrated  in  Figure  tf.  The  detector  noise  is 
assumed  to  be  wideband  in  order  to  model  the  phase  variance  of  the  PLL 
output,  and  is  modeled  as  shown  in  Figure  9.  The  important  result  from 
Figure  9 is  that  the  phase  variance  is  equal  to  1/SNR,  where  SNR  is  the 
SNR  in  the  closed  loop  of  the  PLL,  B (Ref  ' : 1 . (This  is  not  to  ue 

Li 

confused  with  the  bandwidth  of  the  spatial  filter.)  The  phase  variance 
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Fig.  8.  Range  Distortion  Processes  for  the  Laser  Line 
Scanning  Sensor 
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Fig.  9.  Phase  Variance  in  the  PLL 


can  be  directly  related  to  the  range  variance  by  using  the  relation, 
& = ^Trrumr  , Thus,  the  range  variance  can  be  expressed  as  follows: 


'I  - ^ 2 <ISr> 


(Ref  3:15-16)  (IT) 


For  the  model  of  the  phase  variance  of  the  PLL  in  the  receiver,  it  is 
assumed  that  the  PLL  is  operating  in  its  linear  region  which  implies 
that  the  phase  variance  must  be  below  a certain  threshold.  Typically , 
this  threshold  is  chosen  to  be  .25  (Ref  3:16-17),  or  that  the  closed 
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loop  SNR  of  the  detected  signal  must  be  greater  than  four.  From  the 
model  of  the  range  error  caused  by  the  receiver  noise  in  the  PLL-,  it  is 
apparent  from  the  reciprocal  relation  in  Figure  9 that  as  the  SNR 
increases,  the  range  error  decreases.  From  Eq  (IT),  it  is  apparent  that 
as  the  modulation  frequency  of  the  laser  is  increased,  the  range  error 
decreases . 

Using  the  constant  reflectivity  assumption  and  setting  r(x)  equal 
to  a sinusoid,  the  spatial  filter  model  can  be  represented  as  in  Figure 
10.  If  the  Bessel  Series  representation  of  Eq  (IS)  is  used  and  b(x)  is 


Fig.  10.  Single  Sinusoid  Range  Input 
approximat  e i as  an  ideal  low  pass  filter,  then  the  output  can  be  easily 
represented  in  terms  of  the  Bessel  Series  input  as  shown  in  Figure  11, 
where  N equals  the  largest  integer  portion  of  the  bandwidth  of  the  spatial 
filter  divided  by  the  frequency  of  the  sinusoidal  range  variation.  The 
approximation  of  b(x)  as  an  ideal  low  pass  filter  is  a reasonable 
simpli fication  since  the  beam  intensity  function  is  typically  Uaussian 

shaped  which  is  comparable  to  the  function  up  to  its  first  oero 

nl 

The  output  of  the  spatial  filter  is  a finite  Bessel  Series  which  is 
easily  separated  into  real  and  imaginary  parts  by  using  the  Bessel 
Series  representation  of  Eq  (io).  The  output  phase  can  be  expressed  in 


Fig.  11.  Output  Representation  of  the  Spatial  Filter 


the  following  form: 


L^-l^  imaginary  teruis  of  Bessel  Series  [ Eq  (lo)h 
* 'real  terms  of  Bessel  Series  [£4  (lo)]  ' 

This  expression  for  the  output  piiase  of  the  spatial  filter  is  used  to 

compute  the  phase  error  by  subtracting  the  output  pnase  from  the  input 

phase;  however,  the  phase  error  changes  according  to  the  value  of  x.  A 

measure  of  the  phase  distortion  which  accounts  for  these  changes  is  the 

integrated  mean  squared  phase  error  which  is  the  phase  error  squared  and 

averaged  over  the  period  of  the  sinusoid  range  input,  zero  to  fir.  Eq  (19) 

is  the  expression  for  the  integrated  mean  squared  phase  error.  The 

computed  phase  errors  are  limited  to  Mod  to  correspond  to  the  relative 

pnase  measurement  by  the  sensor.  The  mean  squared  range  error  is  computed 


V2  ~ £ /r($(x)  - 0(x))2dx 
l o 
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squared  range  error.  This  is  a convenient  form  for  analysing  the 
distortion  because  it  incorporates  the  maximum  amplitude  variations  of 
the  range.  From  this  discussion  of  spatial  filtering,  it  should  be  noted 
that  as  the  frequency  of  modulation  increases,  so  does  the  range  distor- 
tion error  because  the  phase  error  increases  when  the  argument  of  the 
Bessel  Functions  in  Eq  (18)  increase. 

The  total  mean  squared  error  is  computed  by  adding  the  distortion 
from  the  PLL  and  the  spatial  filtering.  For  low  frequencies  of  modulation, 
the  total  distortion  error  is  dominated  by  the  receiver  noise  distortion 
process,  while  at  high  frequencies  of  inoculation  the  total  distortion  is 
dominated  by  the  spatial  filter  distortion  process.  Somewhere  in  between, 
there  is  a frequency  of  modulation  where  the  range  distortion  is  minimal, 
i.e.,  optimum  modulation  frequency. 

One  other  point  which  needs  to  be  addressed  is  that  the  SNR  depends 
on  the  bandwidth  of  the  PLL.  For  comparative  distortion  calculations , it 
is  assumed  that  the  closed  loop  bandwidth  of  the  FLL  equals  the  bandwidth 
of  the  spatial  filter.  Thereby,  the  spatial  filter  bandwidth  varies 
according  to  the  beam  spot  sine,  and  it  is  necessary  to  normalize  the 
PLL  range  variance  which  is  given  in  Eq  (17)  by  multiplying  it  by  N 
(N  = l/(Xf^)),  where  different  N's  relate  to  different  beam  spot  sizes. 

This  normalized  SNR  will  be  designated  by  CNR 

norm 

Distort  ion  Computations  and  Res  ul t s 

Distortion  results  were  evaluated  by  using  the  sinusoidal  range 
model.  Actual  calculations  were  computed  on  a CYBER  7^  computer,  and 
tiie  Bessel  Functions  were  calculated  by  an  IBM  InO  subroutine  called 
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BESJ.  A brief  list  of  the  assumptions  used  for  making  the  distortion 
computations  are  as  follows: 

1.  The  temporal  scanner  modulation  is  periodic  and  only  a single 
harmonic  is  processed  by  the  receiver. 

2.  The  terrain  reflectivity  is  constant  over  several  beamwidths  or 
its  dependency  has  been  removes  from  the  input  signal  of  the  spatial 

filter. 

3.  The  terrain  range  variation  is  modeled  as  a single  sinusoiu. 

4.  The  beam  function,  b(x),  is  modeled  as  an  ideal  low  pass  filter. 

5-  The  detector  noise  is  wideband. 

6.  The  PLL  is  operating  in  the  linear  region. 

7.  The  loop  bandwidth  of  the  PLL  is  equal  to  the  spatial  filter 

bandwidth. 

The  numberical  integration  of  Eq  (19)  is  implemented  by  using  2000  inter- 
vals of  integration.  For  convenience,  the  distortion  data  is  presented 
re" 

as  plots  of  — r versus  B (the  phase  modulation  index),  where  B is  equal 

* m 

( 7Tn~mrm).  in  all  cases  considered,  only  the  first  harmonic  of  the 
c 

periodic  signal  is  processed  by  the  receiver  (i.e.,  n = l)  ; and  henceforth, 
n will  be  dropped  from  the  expressions  of  B-  Figures  12-18  are  the  plots 

re‘- 

of  — — versus  B for  comparisons  of  different  beam  spot  sizes  (N  as 

rnT 

related  to  f and  X)  and  for  various  SNR  's  and  will  be  discusses 
r norm 

later.  The  absolute  SNR  is  related  to  the  normalized  SNR  by  the  following 

expression : 

SNR  v = SNR  - 10  log  N (20) 

abs  no  rm  ° 

wnere  SNR  , is  the  absolute  SNR.  The  clots  of  the  range  distortion  data 
abs 

can  be  applied  to  distortion  calculations  when  the  loop  bandwidth  of  the 
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RANGE  ERROR  SQUAREO  OIVIOEO  BY  R-MAX  SQUARED 


RANGE  ERROR  SQUARED  01 


PLL  exceeds  the  bandwidth  of  the  spatial  filter  by  using  the  expression 
in  Eq  (20).  These  distortion  results  can  not  be  applied  if  the  loop 
bandwidth  is  less  th:ui  the  spatial  filter  bandwidth  as  there  will  be 


additional  distortion  caused  by  the  PLL  which  acts  as  a low  pass  filter. 

The  curves  show  plots  of  the  normalized  mean-squared  range  distortion 

error  versus  3 for  N = 1,  3,  5,  7»  and  9 • The  even  values  of  N have  been 

deleted  from  the  plots  so  that  the  curves  do  not  become  cluttered, 

however,  the  results  fall  in  between  the  odd  values  of  N as  expected. 

For  values  of  N less  than  one,  the  output  phase  is  zero  and,  thus,  the 

sensor  cannot  resolve  ;uiy  objects  smaller  than  one  beamwidth.  Also,  the 

value  of  N is  always  an  integer,  and  if  the  computed  value  of  N (l/(Xf^)) 

is  not  an  integer,  the  fractional  portion  is  truncated.  The  reader 

should  be  reminded  that  the  data  in  the  plots  is  in  terms  of  relative 

range  error,  not  absolute  range  error. 

The  distortion  data  which  is  plotted  for  various  M's  and  SNR  ' s 

norm 

shows  some  distinct  trends  between  the  sensor  and  terrain  parameters  and 

the  range  distortion  error.  In  analyzing  the  distortion  data  r will 

always  be  fixed  and  ch:uiges  in  R will  represent  changes  in  the  frequency 

of  modulation,  f . In  all  these  curves,  tile  tail  of  the  curves  drops 
ra 

off  for  large  R and  small  N's  because  the  l/r~  normalising  term  is 


dominant , not  because  the  absolute  range  distortion  error  is  getting 
smaller.  As  illustrated  in  tile  tail  of  these  curves , the  range  distor- 
tion error  does  saturate  because  the  phase  error  is  Limited  to  Mod  .n. 

As  previously  predicted,  there  is  an  optimum  operating  value  of  R,  which 
is  interpreted  as  the  opt  Lmum  fro  pienoy  of  modulat  ion  whor  ■ t he  range  err  v 
distortion  is  minimal.  For  example,  if  the  laser  line  scanning  sensor 
is  used  for  mapping  specific  terrain  'file.;  >r  ' ir  •:  s • • frequency 
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of  modulation  could  be  specifically  chosen  so  that  the  minimum  distortion 
of  the  range  signal  occurs  for  that  target  profile.  The  minimum  distor- 
tion point  decreases  as  the  signal  to  noise  ratio  increases.  The  value 


of  3 for  the  minimum  distortion  point  also  decreases  as  the  SNR  increases. 
For  larger  values  of  N (larger  spatial  bandwidths),  the  minimum  distor- 
tion point  is  at  larger  values  of  3 or  higher  f . Based  on  this  analysis 
of  the  curve  trends,  it  is  desirable  to  manipulate  the  modulation  frequency 
so  that  the  sensor  is  operating  at  the  minimum  distortion  point. 

The  question  arises  as  to  how  these  plots  can  be  used  to  determine 
distortion  for  a specific  set  of  parameters.  B is  determined  by  computing 
, UTTfmr,.. . 

K Si.  ■ ) , where  r is  the  peak  amplitude  height  variation  of  the  terrain 

range  profile,  and  f is  ttie  frequency  of  modulation  of  the  laser.  N is 

calculated  by  dividing  the  spatial  bandwidth  of  the  beam,  1/X,  by  the 

spatial  frequency  of  the  terrain,  f . The  number  of  harmonics  at  the 

output  of  the  spatial  filter,  N,  must  be  an  integer.  If  the  absolute  SNR 

is  measured,  then  Eq  (20)  must  be  used  to  calculate  the  SNR  to  use 

norm 

the  appropriate  figure.  A sample  calculation  follows: 

Sample  Calculation 


Given:  Beam  Spot  Sise 

X = 15  cm 


Frequency  of  Modulation 

f = 100  MHs 
m 


Maximum  Height  Variations  oNR  . = 20db 

, ubs 

r = lm 
m 

N = -j— y (truncated  to  an  integer  value)  = r> 
r‘ 

h t f r 

B = m m = h.i9 


Spatial 
Frequency  of 
the  Range 
Variations 
fr  = 1-25 

cycles/m 


SNR  = SNR  . + 10  Log  N 

norm  abs 


SNR 


norm 


ub 


ob 


To  find 


2 


, Figure  17  is  used  and  the  curve  for  N = ‘3  witli 


l 


r 

e 


r 

m 


= 4.19. 


8 x 10 


-4 


r ^ = 8 x 10  m4- 

e 


.03  m (or  about  3%  distortion  of  the  peak  range  amplitude 
variation) 


The  terrain  parameters  f and  r can  be  interpreted  as  the  spatial 

r m 

frequency  and  the  peak  amplitude  variations  of  the  terrain  range  profile, 

respectively.  For  instance,  if  the  profile  of  an  automobile  was  being 

mapped,  the  peak  amplitude  variation  that  would  be  desirable  to  discern 

would  be  approximately  1.5m  and  the  maximum  spatial  frequency  that  would 

be  desirable  to  discern  would  be  10  cycles/m  or  a range  variation  with  a 

period  of  10  cm.  This  is  the  type  of  interpretation  used  to  describe 

the  particular  target  or  range  profile  parameters. 

As  it  is  desirable  to  operate  with  minimum  distortion,  t no  optimum 

frequency  of  modulation  is  plotted  against  the  normalised  signal  to  noise 

ratio  in  Figure  19  for  r = lm.  This  single  plot  with  one  value  for  r 

m ° r m 

is  sufficient  for  describing  the  relation  of  all  other  r 's  because 

in 

the  curve  scales  for  other  values  of  r bv  Lho  factor  1 r . This  is 

m • m 

apparent  from  the  previous  riots  where  the  optimum  modulation  freeucncv 
is  determined  bv  taking  the  value  of  3 at  the  first  minimum  of  the  curve 
and  dividing  the  value  of  £ bv  ( -.'!/-■?.)  , and  noting  that  f is  sealeu  by 
1/r  . The  upoer  limit  on  the  optimum  modulation  frequency  of  >:  10' 

Hz  in  Figure  lo  is  the  result  of  the  largest  values  o*’  3 which  were  used 
in  the  numerical  calculations , not  tho  distortion  process. 
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OPTIMUM  MODULATION  FREQUENCY(HZ) 

S 6789J01  2 3 4 5 6 78910*  2 3 4 5 6 789 


Some  additional  distortion  calculations  X'or  some  specific  terrain 
and  sensor  parameters  are  presented  in  Table  1.  In  Table  1,  various 
sensor  and  terrain  parameters  have  been  varied,  and  the  root  mean  squared 
range  distortion  error  has  been  computed  from  the  previous  curves.  The 
resulting  range  distortion  error  is  a siceable  percentage  of  r in  all 
cases;  thus,  the  distortion  process  severely  limits  the  range  accuracy 
of  the  sensor  for  these  particular  cases.  In  Table  1,  if  case  i is 
compared  to  case  2,  the  spatial  frequency  is  smaller  for  case  2 and  the 
resulting  rms  range  distortion  error  is  smaller.  If  case  1 is  compared 
to  case  3,  where  r is  Larger  in  case  3,  the  rms  range  distortion  error 

Cl 

is  ten  times  larger  in  case  3.  A comparison  of  case  -t  and  case  5 
illustrates  the  effect  of  decreasing  the  beam  spot  sice.  Cases  o-lO  are 
some  examples  where  several  of  the  terrain  and  sensor  parameters  have 
been  varied  at  the  same  time.  Table  2 is  a presentation  of  the  same 
cases  used  in  Table  1,  but  the  optimum  frequency  of  modulation  has  been 
picked  for  the  sensor  by  using  the  curves  presented  in  Figure  Id.  Sy 
comparing  the  root  mean  squared  range  error  in  Table  1 and  Table  2,  for 
most  of  the  eases  there  are  significant  decreases  in  the  range  distortion 
error  when  the  optimum  modulation  frequency  is  used.  Even  so,  the  range 
distortion  still  exceeds  10J*  of  r in  cases  i,  3,  4,  6,  7,  8,  and  10, 

til 

which  means  that  a smaller  beam  spot  sice  will  be  tiie  only  way  to  reduce 
the  distortion  any  further. 

If  the  sensor  is  used  in  a forward  looking  configuration  ;.s  shown  in 
Figure  20,  then  the  distortion  results  can  be  extended  to  evaluate  its 
range  distortion  performance  when  the  system  is  used  at  low  gracing 
angles.  Both  tiie  beamwidth  and  the  range  ;unplitude  variations 
seen  by  the  sensor  are  modified  from  their  original  definitions 


wm  on  are 


Table  2.  Some  Range  Distortion  Calculations  with 
Minimum  Distortion 


Fig.  20.  Forward  Looking  Configuration  for  the  Laser 
Terrain  Mapping  Sensor 


(2l)  the  expressions  for  the  new  r and  X are  presented. 

m 

X'  = X/sin  a 


(21a) 


r ' = (X'  )*•  - (r  sina) 

m m 


(21b) 


This  is  not  a complete  description  of  the  forward  looking  configuration 
as  the  fii'  ■ function  would  also  have  to  be  modified,  but  it  is  reason- 
able for  obtaining  results  for  design  calculations  to  check  the  feasibility. 
A sample  calculation  for  the  forward  looking  configuration  follows: 

Sample  Calculation  for  Forward  Looking  Configuration 


Given:  X = .15  m 

f = 1 x 10"'1  Hz 
m 

r = .2u  m 

m 

■ Calculated:  X'  = .5t)  m 


f =1.25  cycles /m 
r 

SNR  , = SOdb 

abs 

01  = 1*5° 


r'  = .56  m 

m 


r - 5 1.9  x 10  and  ,2  = .025  m 

e r 

r S e 


t*  = S.4i 

ib.5db  (Need  to  interpolate 
between  Figures  1.' 
and  18) 


* 

IV.  Conclusions  and  hoc'’:j’-g:sat  i 

The  basis  ct  this  theoretical  study  was  to  analyze  the  spatial 
filtering  distortion  effects  on  the  performance  of  the  laser  line  scanning 
sensor.  Some  specific  distortion  results  were  presented  wnieh  showed  the 
interaction  between  the  terrain  characteristics  and  the  sensor  parameters. 

From  this  analysis  and  the  distortion  results,  several  conclusions  and 
reoommenuations  can  be  made. 


<,CTlC. Lao  1 sT  l i 


This  study  began  with  an  examination  of  a one  dimensional  spatial 
filtering  model.  An  attempt  was  made  to  model  trie  terrain  range  and 
reflectivity  as  Gaussian  Random  Processes.  This  mouel  of  the  terrain 
only  served  the  spatial  filtering  analysis  by  describing  the  conuitions 
for  no  distortion.  Next,  a detailed  examination  of  the  distortion 
enacts  was  maue  to  determine  what  simplirications  were  necessary  to  maxe 
the  analysis  more  tractable.  Specifically , two  approaches  were  modeled, 
linearization  of  the  phase  modulated  range  signal  ana  the  case  when  the 
terrain  height  profile  was  constant.  In  addition,  a model  for  sampling 
receivers  as  applied  to  pulse  types  of  modulation  for  the  laser  was 


analyzed. 


The  next  ster  in  the  analysis  was  to  obtain  some  distortion  results 
which  were  useful  in  ueterminir.g  the  interaction  between  tiie  sensor  and 
terrain  parameters  :md  their  effect  on  the  distortion.  First,  tiie  output 
of  tiie  spatial  filter  was  represented  in  terms  of  tne  input  terrain 
z urometers  * Anulyticuliy , tais  soxiit  ion  was  sot  very  tractable  and  once 
again  s imp  1 i 1 1 cat  i ons  vein  iv*  u \ . . 1 si  nip  ..  1 1 icat  ion  i u vv  1 vou  the 

removax  of  the  reflectivity  so.  s ieuco  which  reduced  tiie  comoloxity  of 
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the  analysis.  But  the  rrinarv  approach  which  was  ensued,  was  to  model 
the  terrain  height  profile  as  a single  sinusoiu  and  to  assume  that  the 
reflectivity  was  constant  so  that  the  received  signal  was  the  phase 
modulated  terrain  height  profile.  Then,  this  phase  modulated  signal  was 
filtered  by  an  ideal  low  pass  filter  which  represented  the  effects  of 
the  spatial  filter  (finite  beam  sice).  In  addition,  another  distortion 
process,  caused  by  the  receiver  noise  in  the  pnase  tracking  circuitry, 
was  modeled,  using  the  two  distortion  processes,  the  total  range 
distortion  error  was  calculated  for  various  terrain  anu  sensor  parameters, 
and  was  presentee  as  set  of  plots.  Some  specific  distortion  calculations 
ware  mace.  From  the  range  distortion  error  curves  and  calculations, 
several  trends  in  the  data  were  observed,  liven  the  maximum  height 
variation  of  the  terrain,  there  was  an  optimum  frequency  of  modulation 
for  the  laser  where  the  range  distortion  error  was  minimised.  This 
minimum  distortion  point  was  found  to  be  dependent  upon  the  SNK,  the 
spatial  frequency  of  the  terrain,  the  beam  spot  sice,  and  the  maximum 
height  variations  of  the  terrain.  The  optimum  frequency  of  modulation 
decreased  when  the  SNR  increased,  and  the  range  distortion  error  at  this 
point  decreased.  In  addition,  when  the  beam  spot  sice  was  reduced  or  tae 
bandwidth  of  the  spatial  filter  was  increased,  the  optimum  frequency  of 
modulation  decreased,  and  the  range  distortion  error  at  tnis  point 
decreased.  The  range  distortion  error  increased  when  the  spatial 
frequency  of  the  terrain  range  variations  was  increased.  Based  on  these 
trends  and  the  distortion  data,  it  was  determined  that  the  sensor  perfor- 
mance could  be  improved  by  optimising  the  frequency  of  moduiat Lon  of  the 
laser  for  a particular  set  of  terrain  characteristics  ar.d  sensor  para- 
meters. rite  data  was  restructured  as  a clot  of  tae  frequency  of  modula- 


te 


tion  of  the  laser  for  minimum  distortior 


i versus  the  signal-to-noise 
ratio.  Finally,  on  extension  of  the  distortion  calculations  was  applied 
to  the  forward  looking  configuration  of  the  terrain  mapping  laser  sensor. 

The  following  conclusions  can  be  drawn  from  this  thesis: 

1.  The  performance  of  the  laser  line  scanning  sensor  can  be 
improved  by  removing  the  reflectivity  dependence. 

2.  For  pulsed  types  of  scanner  modulation,  it  was  determined  that 
the  output  of  the  spatial  filter  was  bandlinited  by  the  bandwidth  of  the 
spatial  filter. 

3.  Throughout  the  analysis  no  information  was  available  about 
terrain  features  which  were  smaller  than  the  beam  spot  size. 

d.  The  performance  of  the  terrain  mapping  sensor  could  be  improved 
by  selecting  the  optimum  modulation  frequency  for  a given  set  of  terrain 
and  sensor  parameters. 

Recommendations 

Several  recommendations  resulted  from  this  study  and  all  of  them 
apply  to  follow-on  studies  of  the  sp>atial  filtering  distortion  problem. 
First,  the  spatial  filtering  distortion  should  be  re-examined  and  recalcu- 
lated for  other  periodic  and  nonperiodic  range  functions  using  other 
distortion  measures  and  calculation  techniques  from  different  sources 
(which  are  included  in  the  Bibliography).  The  distortion  range  error 
should  be  recalculated  for  a Gaussian  shaped  filter  which  more  closely 
represents  the  spatial  filter.  Furthermore,  the  results  which  were  ob- 
' tair.ed  do  not  apply  to  time  domain  processing  receivers  where  all  the 
harmonics  of  the  periodic  modulation  are  detected.  For  these  calculations. 
Fast  Fourier  Transform  techniques  will  probably  serve  as  a useful  numerical 
technique,  and  eouid  be  utilized  tc  analyze  specific  height  profiles  of 
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